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ABSTRACT 

Observations of higher-excited transitions of abundant molecules such as CO are important for determining where energy in the 
form of shocks is fed back into the parental envelope of forming stars. The nearby prototypical and protobinary low-mass hot core, 
IRAS 16293-2422 (116293) is ideal for such a study. The source was targeted with ALMA for science verification purposes in band 
9, which includes CO J=6-5 (£^upAb~116 K), at an unprecedented spatial resolution (~0V2, 25 AU). 116293 itself is composed of 
two sources, A and B, with a projected distance of 5''. CO J=6-5 emission is detected throughout the region, particularly in small, 
arcsecond- sized hotspots, where the outflow interacts with the envelope. The observations only recover a fraction of the emission in 
the line wings when compared to data from single-dish telescopes, with a higher fraction of emission recovered at higher velocities. 
The very high angular resolution of these new data reveal that a bow shock from source A coincides, in the plane of the sky, with 
the position of source B. Source B, on the other hand, does not show current outflow activity. In this region, outflow entrainment 
takes place over large spatial scales, >100 AU, and in small discrete knots. This unique dataset shows that the combination of a 
high-temperature tracer (e.g., CO J=6-5) and very high angular resolution observations is crucial for interpreting the structure of the 
warm inner environment of low-mass protostars. 
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IRAS 16293-2422 (hereafter 116293) is a well-studied solar- 
type protostar, originally singled out as having one of the most 
prominent outflows of all the IRAS sources associat ed with the p 
Ophiuchus cloud, and showing indications of infall ( [Walker etal.| 
|1986| |1988| ). Sub-millimetre and radio observations reveal that 
the source breaks up into two sources, A and B, lo cated at a pro- 
jected distance of 5'' ( [Mundy et al.|1986[|Wootten| 1989 ). Source 
A itself breaks up into two sources, Al and A2 separated by 0V3 
([Wootten|1989| ). At the distance to the source (120 pc, [Loinard| 



et al. 2008), these angular separations correspond to physical 
distances of -600 AU and 36 AU respectively. While sources 
Al and A2 appear active and drive two outflows, source B is 



I more quiescent. Infall is clearly detected towards source B (e.g. 
[Pineda et al.|2012|) and several sho ck tracers such as SiO are de- 



tected (e.g., J0rg ensen et al.|2011[ ), a lthough other obs ervations 
indicate source B is a prestella r core ([Chandler et aTp005 ) or a 
T Tauri star ( [Stark et al.|2Q04l[Takakuwa et al.|2007t . Figure p] 



shows the locations of the diflTerent protostellar components 

Of the two sources A and B, source B has been the more puz- 
zling because molecular lines observed towards the so urce show 



a large spread in Av ranging from ~ 1 to 5 km s (J0rgensen 



[et al.|20lT] ). Towards source A, on the other hand, the lines all 
have a narrow range of Av of ~ 2-3 km s"^ A high angular 
resolution molecular line survey carried out with the Submil- 
limeter Array (SMA) shows that several species peak towards 
both sources whereas some are only detected towards source A, 



* Present address: Harvard-Smithsonian Center for Astrophysics, 
60 Garden Street, MS78, Cambridge, MA 02138, USA, e-mail: 
lkristensen@cf a . harvard . edu 



others only towards B ( jBisschop et al. [[20081 [J0rgensen et al. 
201 1[ ). ALMA observations towards source B have recently 
shown infall in CH3OCHO ( [Pineda et al.|2012[ ) and the first de- 
tection of the simplest sugar, glycolaldehyde, as well as several 
related molecules (J0rgense n et al.[[2012[ ). The diflTerent mor- 
phologies and velocity structures revealed by diflTerent species 
make 116293 one of the most enigmatic low-mass protostellar 
sources as illustrated by the wealth of data available. 

Understanding a protostellar system such as 116293 requires 
not only an understanding of the inf ailing envelopes, but also 
when and where energy is released back into the envelope, pri- 
marily in the form of shocks driven by the protostellar outflow 
and jet. On large scales (> 10^ AU) 116293 drives a quadrupolar 
outflow ( [Mizuno et al.|1990[ [Stark et al.|2004"| ), but disentangling 
the eflTects of the driving jet on small scales is not straightfor- 
ward. Shock hotspots created by the outflow have the potential to 
change the local physical conditions (density and temperature), 
and are best traced by species which trace currently shocked gas 
and are unaflTected by chemistry. 

In the following letter, we present ALMA science verifica- 
tion data of the CO J=6-5 transition towards 116293. The very 
high angular resolution of these Band 9 observations (0V2 at 
690 GHz), and a high-temperature tracer only weakly aflfected 
by chemistry are a perfect combination to re-examine the origin 
of molecular emission towards these two sources. Although CO 
J=6-5 also traces swept-up outflow gas on larger spatial scales 
(> 1000 AU; [van Kempen"eraLl[2QQ9l [Yildiz et al.[[20T? ), this 
emission is filtered out by the interferometer and therefore these 
new data show only the small-scale structure where the outflow 
is interacting with the envelope. 
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Fig. 1. Left: Map of CO J=6-5 emission towards 116293 with contours at 6, 12, 18, .cr of integrated intensity. Integration limits are -10 
to +lkm s~^ and +7 to +18 km s~^ for the blue- and red-shifted emission, respectively. The underlying gray-scale image shows continuum 
emission at 690 GHz. The positions of sources A and B are marked; source A consists of the two sources Al and A2. Offsets are recorded from 
16^32"'22?753; -24°28'34r747 (J2000). The ALMA beam is shown in the lower right corner (black). The centre of the JCMT beam is marked 
with a cross and the beam is shown as a dashed circle. The dashed lines on source A shows the region used to construct the pv diagram (Fig. |4|. 
Right: The emission integrated over the highest velocities (from -10 to -4 km s~^ and +12 to +18 km s~^) is shown in contours at 5, 10, 15, . . . cr. 
The different features are highlighted by arrows and are labeled. The black dashed arrow is an extrapolation of the red lobe of the NW-SE outflow. 



2. Observations and results 

The observations consist of a seven-point mosaic in Band 9, and 
were taken as one of the ALMA Science Verification datasets 
on April 16 (1 execution) and 17, 2012 (3 executions) using 
15 antennas and a configuration with baselines ranging from 26 
to 402 m (62 to 943 k.1 at 690 GHz). The on-source integration 
time was ~18 minutes per pointing and the primary beam is 9''. 
All four spectral windows consist of 3840 channels with a chan- 
nel separation of 488 kHz (0.21 km s"^ at 690 GHz) and a total 
bandwidth of 1.875 GHz. The calibrated data were taken from 
the CASA guides websit^ We performed phase and amplitude 
self calibration on the continuum, confirming the results from 
the CASA guide. The gain solution from the self calibration 
was then applied to the measurement set, and the spectral region 
around the CO 6-5 line at 691.473 GHz was extracted. 

Due to the combination of large-scale emission of CO and 
the missing short spacings, special care was taken to identify 
and flag data points with amplitude outliers. These could result 
in image artifacts being interpreted as real structures due to in- 
complete sampling of the wv-plane. Since the wv-coverage of the 
dataset obtained on April 16 has a large overlap with one of the 
datasets from April 17, the derived amplitudes were compared 
to help identify bad data. Two out of 105 baselines exhibited 
questionable amplitudes and were therefore flagged. The details 
of this procedure are discussed in Appendix |A| The visibilities 
were imaged using a robust weighting of 0.5. The resulting spa- 
tial resolution is 0V29x0V16 at a position angle of 110°. The 
rms noise is lOOmJy beam"^ per 0.21 km s"^ channel. The spec- 
tral cube and the continuum image were corrected for primary- 
beam attenuation, and smoothed to an effective resolution of 
0.84 kms-^ 

CO J=6-5 emission is clearly detected and the integrated 
emission shows complex structure. Figure [T] shows CO J=6-5 
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Fig. 2. ALMA (black) and JCMT (blue) continuum- subtracted CO 
J=6-5 spectra compared in the same 8'' beam marked in Fig. [T] The 
dashed lines indicate the velocity interval excluded in this analysis 
(i^source ±3 kms~^) and the dotted lines illustrate the high- velocity inter- 
vals less affected by filtering. The horizontal line (green) is the baseline. 



emission overlaid as contours on the continuum, where the blue 
and red- shifted emission is integrated separately from -10 to 
-hi km s"^ and +7 to 18km s"^ respectively (Fig. [2]). As can 
be seen from Fig. [T] there is very little overlap between blue- 
and red-shifted emission except around source A where there is 
some overlap. 

To quantify the amount of spatial filtering in these data, we 
compare emission to a CO J=6-5 spectrum from the James Clerk 
Maxwell Telescope (JCMT) towards I16293-A. The spectrum 
was obtained under good atmospheric conditions in July 2000 
and is similar to that presented in |Stark et al.| ( [2004"| ). The JCMT 
was pointed to a position approximately V N of source A and 
the beam- size is 8'' (Fig. [T]). The JCMT data were converted 
to flux density using a conversion factor of 34.0 Jybeam"^ K"^ 
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Fig. 3. Continuum- subtracted spectra towards the sources A and B. 
The spectrum towards A has been shifted by 2 Jy beam"^ Basehnes 
are shown as full lines, and the velocity range excluded from analysis 
is shown with vertical dashed lines. Unidentified features are marked 
with a U (see text). 



The ALMA data were convolved with an 8'' beam centred at the 
same position (Fig. [2]). ALMA recovers ~5% of the flux density 
in the full velocity range, but when considering only the blue or 
red wings (excluding the central 6 km s"^ from f;LSR ± 3 km s"^), 
approximately 15 and 10% of emission is recovered, respec- 
tively. If the velocity cut-ofl' is moved from f^source ±3 km s"^ 
to ± 8 km s"\ ~ 50 and 15% of emission is recovered in the blue 
and red wings, respectively. These high-velocity structures are 
therefore less afl'ected by filtering. They are shown in Fig. [T] 
(right) where they likely form the edges of larger- scale structures 
not recovered. 

The peak flux densities are typically of the order 1-2 Jy per 
ALMA beam. Figure [3] shows single-point spectra towards the 
continuum peaks of sources A and B. Towards source A, the 
profile is dominated by blue-shifted emission and absorption at 
the source velocity. Towards source B, the CO J=6-5 profile is 
seen almost entirely in absorption against the strong continuum. 
More features are seen in absorption towards source B, similar to 
the observations of glycolaldehyde ( J0rgensen et al.|[2012j ; we 
do not speculate on the nature of the species responsible for the 
two absorption features labeled "U" in Fig. |3] 

3. Discussion: origin of the emission 

The highest- velocity emission is shown in the right panel of 
Fig. [T] where it outlines active hotspots where the outflow is in- 
teracting with the envelope. We therefore focus the discussion 
on their morphology and the interpretation. Source A hosts a 
red-shifted arch capped by blue-shifted emission pointing NW 
and blue-shifted emission surrounded by red-shifted emission 
stretching SE. The wing emission is strongest and most extended 
in terms of velocity towards the region SE of source A. Near 
source B an arch of blue-shifted emission opens to the SE to- 
wards source A. 

Concerning the knots around source A, the origin of the 
large-scale E-W outflow is seen in two outflow knots that are 
aligned with source A2. This is consistent with the observations 



of |Pech et aL| ( |2QT0| ) and |Loinard eTaLl ( [20T2l ). On these smaU 
scales the outflow direction is switched in the sense that the red- 
shifted lobe is moving west, while the blue-shifted is moving 
east, oppos ite to what is observed on large scales (e.g., [Stark| k 
et al.|2004] ). The lobes of the large-scale E-W flow do not extend " 



continuously to the inner region, implying that the system is en 
tering a new accretio n phase with a slightly diflTerent inclination 
angle ( |Pech et al.|2010 , Loinard et al.|2012| ) . 

The distance traversed by the red- shifted shock of the E-W 
flow (the tip of the eastern red arrow in Fig. [T]), ~0V5, since 
its launch in 2005 is consistent with proper motions seen in the 
cm ob servations of this outflow knot (~0V1 yr~^; |Pech et al. 
2010) and implies a space velocity of ~90km s"^ with an in- ^ 
clination angle of -5° with respect to the plane of the sky. With 
ALMA it is possible to trace the motion of these outflow knots 
on timescales of a few years. Thus, the velocity and momen- 
tum transfer can be followed in 3D, providin g tight constraints 
on entrainment simulations ( Arce et al.|2007| ). 

There is a clear asymmetry between the blue- and red- shifted I 
emission of the NW-SE outflow in terms of spatial distribution. 
The red-shifted emission NW of source A breaks into 3-4 knots 
that are aligned in a semi-shell. These knots may represent small 
shock condensations where the outflow is directly interacting 
with the envelope, as is also prominently traced by thermal H2O 
emission towards other sources ( Kristensen et al. 2012 ) and H2O 
maser emission towards this source ( Woott en et al . 1999). The 
blue-shifted emission SE of source A however, shows a highly 
collimated structure bracketed by red-shifted emission. 

To examine whether the cause of the asymmetry is internal 
(e.g., asymmetric launch) or external (e.g., density gradient or 
clumps) we created a position-velocity (pv) diagram with a wide 
slit of 1V25 to capture the red-shifted emission NW of source 
A (Fig. [4]). The line profiles extend out to ±20 km s"^ in these 
knots. The velocity of the blue-shifted emission is highest at the 

base of the outflow (v 20 km s"^), closest to the source but 

rapidly decelerates over 0V5 (~ 60 AU) to just a few km s"^ Part 
of the red- shifted emission appears to be hidden by the contin- 
uum and is only in emission at a distance of ~ 0V3 (35 AU). The 
maximum velocity is also 20 km s"^ At a distance of (120 
AU) from the source the velocity of the red-shifted emission has 
also dropped to the ambient velocity. The symmetry in the pv 
diagram suggests that the origin of emission is the same and the 
cause for the asymmetry is in the surrounding environment, but 
on small scales (< 25 AU). 

The rapid deceleration seen in the pv diagram is different to 



the Hubble-type flows observed over larger spatial scales ( |Arce 
|et al.| r2007 , and references ther ein) but more cons istent with a 
"turbulent wind"-type scenario ( [Canto & R aga 1991 ) where the 
outflow material decelerates through interaction with the sur- 
rounding envelope. Furthermore, the pv diagram illustrates that 
entrainment takes place over scales larger than the inner -100 
AU, at least in this system, and in localised hotspots. 

The arch near source B delineates an empty half- shell; the 
arch likely contains large-scale filtered emission — full ALMA 
(with the Atacama Compact Array, ACA, and Total Power, TP, 
telescopes included) will be required to capture this emission. 
This arch is a coherent structure with a characteristic line width 
of ~ 5 km s"^ which is more consistent with the line widths of 
source A than B. The arch is seen in absorption against source 
B (Fig. [3]), which is oflfset from the apex of the arch. The oflf- 
set is particularly clear in individual channel maps (Fig. |B.4| 
in Appendix |B]). The arch is aligned with the NW-SE outflow 
originating in source Al (Fig. [T]). If source B were the driv- 
ing source of this structure we would expect the highest- velocity 1 
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Fig. 4. Position-velocity diagram of a slice running at a position 
angle of 128° and with a width of 1V25 covering the central part of the 
outflow centred on source Al. The velocity has been shifted so the 
source velocity is at km s~^ (black line). Signal below the 3cr level 
is ignored as are the central +3 km s~^ The position of source Al is 
marked with a dotted line. Negative offsets are towards the SE, positive 
towards NW. 



gas to be closest to source B, which is not the case. If source B 
were the origin of the structure and the blue arch were an outflow 
cone, the outflow structure at the base of the outflow within 25 
AU would be significantly diff'erent from that of source A (see 
Appendix [B] for a derivatio n). Contrary to the very recent inter- 
pretation byfLoinard et al.|('2012), these facts are consistent with 
an origin in source A, and that the arch is a bow shock originat- 
ing from jet interacting with the ambient material surrounding 
source B. 

At the base of the outflow near source A itself, both red- 
and blue-shifted emission is observed and overlapping as can 
be seen both in the map (Eig. [T]) and the pv diagram (Eig. |4]), 
indicating that the outflow is moving close to the plane of the 
sky. [Chandler et al .| p005) reached a similar conclusion based 
on the large-scale outflow maps of (Stark et al.| ( |2Q04| . Even 
small-scale precession in and out of the plane of the sky could 
cause the arch to be blue-shifted at a distance of ~ 5'' (600 AU). 
The dynamic age of the blue arch is < 300 yrs for a flow velocity 
of 10 km s-i (Eig. [3]), an upper limit on the age because the 
flow velocity is a lower limit. Therefore, any precession in and 
out of the plane will have to be on timescales shorter than ~ 
300 yrs. | Pech et al.|([20T0|) demonstrate that Al and A2 are a 
binary system and Chandler et al.| ( [2005| find that the position 
angle between the two sources varies by 2.2° yr"^ whereas the 
projected distance remains constant, implying an orbital period 
of ~ 165 yrs. It is therefore plausible that the binary interaction 
is responsible for some of the diff'erence in motion with respect 
to the plane of the sky. 

Several species have been detected inside the arch with the 
Submillimeter Array (SMA) ([Chandler et al.||20Q5l [Takakuwa 
et al.|2OO7|[j0rgensen et al.|2 011), including SiO and SO which 
are both known shock tracers. The SiO channel maps by !J0r-| 
gensen et al.| pOllj ) show that the SiO emission is both blue- 
and red-shifted, and SMA maps of CO J=2-l and J=3-2 show 



that there is significant blue- shifted emission towards source B 
dTakakuwa et al.[ |20071 [Yeh et al.|[20Q8| ). These observations 
corroborate the scenario that the arch is an active shock region. 
This interpretation is also consistent with the large range of Av 
observed towards source B: some of the species are not asso- 
ciated with this source but arise in the shock originating from 
source A along the same sightline. In particular, the lines with a 
width of ~1 km s"^ likely originate in the inner protostellar en- 
velope, whereas the broader lines (~5km s"^) originate in this 
bow shock. Source B itself appears to be quiescent at the mo- 
ment and is not actively accreting although larger- scale infall is 
observed towards the source ( [Pineda et al.|2012| ). 

4. Summary and conclusions 

CO J=6-5 is detected in a large number of distinct clumps to- 
wards 116293 where shocks are likely active. More flux density 
is recovered in the high-velocity line wings by the interferome- 
ter, and so it is not surprising that only localised higher- velocity 
hotspots are detected. The new data show that there is no direct 
evidence for current shock activity associated with source B, and 
that although infall is observed towards the source, it is not cur- 
rently accreting. Instead, a blue- shifted bow shock from source 
A is overlapping with source B in the plane of the sky. Outflow 
entrainment takes place over large scales, >100 AU, and wind 
material is decelerated through direct interaction with the enve- 
lope. These ALMA science verification data already reveal the 
potential of this new facility for shedding light on how protostars 
interact with the dense inner envelope on spatial scales of < 100 
AU and where energy is fed back into the envelope. 
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Fig. A.l. Channel maps at -3 km s~^ before 
and after additional flagging (see text). The off- 
sets are as in Fig. [T] The color scale ranges 
from to 0.3 Jy beam~^ (0 to ~ 3cr) and has 
been chosen to highlight the striping present 
in the left map (original data) but not the right 
(flagged data). The blue contours are at 1 cr in 
this 0.25 km s~^ channel. The red dot at (0V75, 
1V65 ) sho ws the position where the spectrum in 
Fig. |A.2| was extracted. Sources A and B are 
marked in green. 



Appendix A: Data reduction details 

For extended emission — as seen in CO J=6-5 towards 
IRAS 16293 — the uv sampling density towards short baselines 
is crucial for a successful recovery of large spatial scales. In the 
image reconstruction process, each baseline contributes a corru- 
gation pattern along a position angle in the image plane that cor- 
responds to the position angle of the antenna pair in the uv plane. 
Amplitude and position of this wave-like pattern are determined 
by its complex visibility. A single baseline, therefore, does not 
provide enough information to constrain the spatial origin of the 
emission. The combined information of many antenna pairs with 
various baseline lengths and position angles is required to get an 
image with genuine structures. 

This also means that if a certain region in the uv plane is 
under- sampled, this will lead to imaging artifacts. To a certain 
degree these can be removed during the image deconvolution 
process. However, if an under-sampled baseline exhibits strong 
correlated flux density, the corrugation pattern, which is caused 
by this antenna pair, will remain in the deconvolved image and 
might be misinterpreted as real structure. For this reason we flag 
the baseline DA47&DV14, which is the only one that samples 
spatial scales above 3'\ Since this baseline exhibits high corre- 
lated flux densities, the reconstructed image shows a clear cor- 
rugation pattern, which is indicative of extended structures that 



cannot be quantified correctly (Fig. A.l 



Additional analysis of the visibilities lead to the flagging of 
the antenna pair DA41&DV09. We made use of the fact that two 
execution blocks were observed almost exactly 24 hours apart, 
leading to almost identical uv coverage. This made it possible 
to cross check the visibilities. In certain channels at blue- shifted 
velocities, DA41&DV09 exhibited strong correlated flux densi- 
ties, but with significant diff'erences between the two execution 
blocks. 

The flagging of these two baselines significantly changed the 
spectra at certain positions in the sky (see Fig. |A.2| for an exam- 
ple), particularly at positions associated with the blue arch. We 
suspect that large-scale structure (>3'0 is associated with the 
blue arch, but that full ALMA is required to recover it (with the 
Atacama Compact Array, ACA, and Total Power, TP, telescopes 
included). We note that the overall morphology of the system has 
not changed, and all key structures can be recovered although the 



^ 1.0 




Fig. A.2. Spe ctrum towards a position located along the blue arch and 
marked in Fig. |A.1| The flagged and original data are shown in red and 
blue, respectively. 



blue arch appears to sit on a platea u of e mission at velocities of 
~3 km s"^ in the original data (Fig. |A.l| . 

Appendix B: Comparison of outflow emission at 
sources A and B 

We illustrate below that if the arch were caused by source B, 
instead of being a chance alignment, a very diff'erent outflow 
collimating structure would be required to that found in source 
A. We do so by assuming that the opacity towards the continuum 
peaks of sources A and B are the same (r^ = = r) and that 
the line emission source functions are the same {S^ = = S y) 
towards both continuum peaks, but that only a fraction, /, of the 
emission will be recovered by the interferometer. This fraction 
may be difl'erent towards the two sources, and we denote the 
fraction recovered towards source A with / and towards source B 
with rif. Only in the case where the outflow structure at the base 
of the outflow (i.e., at the continuum peak) is similar towards the 
two sources is = 1 . 
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The standard equation of radiative transfer reads: 
/, =7,(0) e""" +5,(1 -e"^") 



(B.l) 
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where /y(0) is the continuum contribution. Towards sources A 
and B, respectively, the radiative-transfer equations become: 



(B.2) 
(B.3) 



When moving into the blue-shifted Hne wing towards source A, 
the Hne emission decreases until almost nothing is left but the 
continuum: 



(B.4) 



In other words, line emission and absorption are balanced and it 
mH follows from Eq. |B.4| that 

I^(0)(l-e-'^)^fSAl-e-^^). (B.5) 

The ratio of observed emission towards sources B and A is 



7^ ^ lf(0)e--^^r]fSAl-e-'^) 
7A 7A(0)^-v+/^, (1-^-0 

= 7?^"^ + 7/(1 -^-^0 



(B.6) 

(B.7) 
(B.8) 



where R = 7^(0)/7^(0) is the ratio of the continuum emis- 
I sion towards the two peaks. The value of R is 3.05 Jy beam" ^ 
/ 1.13 Jy beam"^ = 2.69. Solving for r gives 



T = In 



R-T] 



(B.9) 



The value of r can be estimated towards source B based on the 
available data. If the absorption is caused entirely by a fore- 
ground layer, r = ln(7v(0)/7v) which is shown in Fig. |B.1[ How- 



ever, the arch is primarily observed in emission around source B 
which means that this value of r is a lower limit because both 
emission and absorption needs to be accounted for. Towards the 
deepest absorption feature, r is greater than 1. 0-1. 5. The ob- 
served intensity ratio, 7^ /I^ is displayed in Fig. 



B.2 



tom of the absorption feature, the emission ratio is 



At the bot- 
0.5. Figure 

B.3| shows 7/ as a function of r for various values of the emis- 
sion ratio 7^/7^. For an observed emission ratio of 0.5, t] is at 
most 0.5, irrespective of r, which implies that at most ~ 50% of 
the emission recovered towards source A is recovered towards 
source B at this velocity. If more emission were recovered, the 
m absorption feature would be filled with the emission. 

The above derivation is only valid if the bases of the outflows 
are similar. As |Loinard et al.| (pOT2|) argue, source B may be 
a candidate first hydrostatic core in which case the outflow is 
expected to be less collimated and therefore may be subject to 
more filtering. However, for that to be the case, source B must 
be located at the apex of the arch. Figure |B.4| shows a set of 
CO channel maps obtained at blue-shifted velocities. All four 
maps, and in particular those at -1.5 and -2.5 km s"\ show that 

• source B is not at the apex of the arch. From these arguments 
we conclude that the blue arch is unrelated to source B and is a 
chance alignment in the plane of the sky. 



i 

I 
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Fig. B.l. Line opacity, r towards the continuum peak of source B. The 
unidentified line features are marked with U. The opacity is a lower 
limit as discussed in the text. 
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Fig. B.2. Line emission ratio from the peak of continuum emission 
towards source B and a line-free position near the peak of source A. The 
unidentified line features are marked with U. 




Fig. B.3. 7/ as a function of the line opacity for diff'erent emission 
ratios. The contours are at emission ratios of 0.2, 0.4, 0.6, ... to 2.0 
with the four lowest ratios labeled. 



L.E. Kristensen et al.: ALMA CO J = 6-5 observations of IRAS 16293-2422 
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Fig. B.4. Channel maps of blue-shifted emission at the velocities 
indicated in the upper left corner of each panel. The channels have a 
width of 0.25 km s~^ and contours are at 4, 8, 12, . . . cr. The background 
grayscale image is the continuum. The positions of sources A and B are 
marked. 
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